Non-Abelian Generalization of Off-Diagonal Geometric Phases 
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If a quantum system evolves in a noncyclic fashion the corresponding geometric phase or holonomy 
may not be fully defined. Off-diagonal geometric phases have been developed to deal with such cases. 
Here, we generalize these phases to the non-Abelian case, by introducing off-diagonal holonomies that 
involve evolution of more than one subspace of the underlying Hilbert space. Physical realizations 
of the off-diagonal holonomies in adiabatic evolution and interferometry are put forward. 
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I. INTRODUCTION 

A quantal system that fails to return to its initial state 
after some prescribed elapse of time may acquire a well- 
defined geometric phase An interesting feature of 
this noncyclic geometric phase is that it becomes unde- 
fined when the initial and final states are orthogonal. 
This gives rise to a nodal point structure that can be 
monitored experimentally in a history-dependent man- 
ner 0, H, 0, [j| . In the hope to recover some of the lost 
interference information at the nodal points of the non- 
cyclic geometric phase, Manini and Pistolesi @ intro- 
duced off-diagonal geometric phases for adiabatic evolu- 
tions of pure states. These quantities may be defined in 
cases where the standard geometric phase is not. The 
adiabatic requirement on the evolution in ref. Q was 
lifted by Mukunda et. al. Q , and Hasegawa et. al. [1, 
provided an experimental verification of the second order 
off-diagonal geometric phase for neutron spin. Theories 
for off-diagonal phases and holonomies for mixed quantal 
states have been developed [H, EH, E2, EH • 

Wilczek and Zee [l4[ showed that the geometric phase 
factor generalizes to a unitary state change, often referred 
to as a non-Abelian quantum holonomy, when consider- 
ing cyclic adiabatic evolution governed by a degenerate 
Hamiltonian. The relevance of non-Abelian holonomies 
for universal fault tolerant quantum computation has 
been demonstrated in refs. [11 El, El El, 13 HI . The 
non-Abelian quantum holonomies have been generalized 
to nonadiabatic [2l|, discrete [221 ]. and noncyclic evolu- 
tions [H, [24J. As for the geometric phase, the holonomy 
may be undefined when the evolution is noncyclic. In the 
non-Abelian case we also have the additional possibility 
that the holonomy is partially defined [24j |. 

In this Letter, we extend ref. [|| and introduce non- 
Abelian off-diagonal holonomies. We demonstrate that 
the off-diagonal holonomies retain holonomy information 



when the standard noncyclic ones [23, [24( are undefined. 
We also suggest physical realizations of the off-diagonal 
holonomies using interferometry, both in adiabatic and 
nonadiabatic settings. 



II. OFF-DIAGONAL HOLONOMIES 



Consider a smoothly parameterized decomposition 



H = n 1 (s)®---®n n (s), se [o,i], 



(i) 



of an A^-dimensional Hilbert space TL into r\ mutually or- 
thogonal subspaces. It is to be noted that this parameter- 
dependent decomposition can arise in different ways, e.g., 
it may consist of the instantaneous eigenspaces of some 
s-dependent Hamiltonian or of some arbitrary de- 
composition that evolves under the Schrodinger equation 
plj . Assume that dim[W/(s)] = n u Vs G [0,1], and 
I = 1, . . . ,r). Thus, each family TLi{s) of subspaces de- 
fines a curve Ci in the Grassmann manifold Q{N; ni), i.e., 
the set of n;-dimensional subspaces in the ./V-dimcnsional 
Hilbert space (2f|. For each such curve, we introduce the 
quantities 



T ; = lim Pi(l)Pi(l - Ss) 



■ Pi(Ss)Pi(0), 



(2) 



where Pi(s) is the projection operator onto the subspace 
H.i(s). This definition makes Ti explicitly gauge invari- 
ant. Note also that the limit Ss — > in Eq. ([2]) makes T; 
uniquely determined for any sufficiently smooth curve Ci 
in the Grassman manifold. We furthermore define 



era =Pk{0)T t 



(3) 
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Let {|/c l (s))}2i and {\l z (s))}Zi be orthonormal bases 
for subspaces 7ik{s) and 7i/(s), respectively, in terms of 
which 

oh = E[^ol^)Pe / ° lAi(s)ds ]y|fe i (0))(F-(0)| 

ij 

= E[^]^l fcl (°))( F (°)l- ( 4 ) 



2 



Here, {J-^\J-[) is a nt x n; matrix with components 
{F&Ttiij = <^(0)|P(1)) and = (W«)|1 J '(*)> 

is the Wilczek-Zee connection along Ci in Q(N; nj). 

The unitary part $[er"] [31] is the holonomy in ref. [2~j ] 
associated with the (open) path C;. It seems natural to 
ask whether we can interpret the matrices tr , k ^ I, 
in a similar fashion. To answer this, we need to see how 
these matrices behave under a gauge transformation, i.e., 
a change of frames \l*{a)) - = l* J '0»)>[tfl(«)U 

where {Ui(s)}^ =1 are unitary matrices. Under such a 

transformation the matrices Pe^o Al ^ ds and (^ol-^i) 
undergo the following changes 

p e tiM s )ds _^ u\(l)Peti Al{s)ds Ui(0), 

(Folrt) - Ulm^iWiH). (5) 
Consequently, cr kl transforms as 

* kl ^U\ (0)o- w E7,(0), (6) 

i.e., noncovariantly unless k — I. Thus, the matrices er fci , 
k I, fail to reflect the geometry of the paths and 
C;. However, the specific behavior of cr' ci under gauge 
transformations suggests that we consider the operator 

7*1. ..2* = a hl K a l K l K ~i ' ' ' a hh a hh 

= E^ 1 "' , "]«l i i(°))< i i( )l' ^ 

where -y' 1 •••'>« is the matrix 

7 /i-2 K = a hl Kcr lJ K -i . .. a hl 2(T hh_ ( 8 ) 

We can use these operators and matrices to de- 
fine gauge covariant quantities, since <i> '"'"] — > 
Lf| i (0)$[7 ili2 - '"][/ ;i (0) from eq. ©. Thus, we propose 
to take 

C/W[C il) ...,C i J = $[V 1 -^] (9) 

as the gauge covariant non-Abelian holonomies of order 
k, and thus generalizing the approach of ref. Q to the 
non-Abelian case. We extend the range of k by defining 
Uy[Ci] = $[er"] , i.e., the first order (k — 1) holonomies 
are taken to be the open path holonomies in refs. [23l |24T| . 

Note that the definition in eq. (J9j) allows any se- 
quence (Zi, •..,/«). This includes cases like, e.g., 7 111 , 
which cannot be regarded as an "off-diagonal" object. 
Hence, eq. (0) can be regarded as a general definition of 
holonomies of degree n, both diagonal and off-diagonal. 
To define genuinely off-diagonal holonomies we obtain 
a reasonable subclass if we require that . . . , l K ) con- 
tains each number at most once. We let 1^ denote 
all vectors (Zi,...,Z K ) with lj G {l,...,r/}, such that 
none of the numbers occurs twice, e.g., (2,5,3) e l| 
but (6,4,6,2) ^ I4. We refer to the set of holonomies 

U^[C h ,...,dJ with (l u ...,l K ) G P K with 2 < K < n, 
as "strictly off-diagonal holonomies" . 



For a cyclic evolution, characterized by = 
H;(0), I = 1,...,T), the standard holonomies U^[C{\ are 
fully defined. On the other hand, in this case we have 
yi...2 K = q 5 ...,/„) g I 1 ^, k > 2, which implies that all 
strictly off-diagonal holonomies are undefined for cyclic 
evolution. Thus, just as in the Abelian case , the stan- 
dard holonomies contain all nontrivial information about 
Ci, . . . , Crf when these are loops. 

In the case where ni = 1, I = 1 , . . . , 77, the matrices 
(J-q \ J-[) and ~PeJo A i( s ) ds reduce to the complex numbers 
(fc(0)|Z(l)) and e -Io(Ks)\d s i(s))ds^ reS p ec tively. This leads 
to the off-diagonal geometric phase factors 

uM[c ll ,...AJ = *[(h(oW)) 

xe-J"o< , -Wl s . , -W>*»... (13(0)11! (1)) 

xe -/o('i(s)l^'i(s))^] j (10) 

which coincide with ^['yl^ 1 ] _ j ] in ref. @. 

Manini and Pistolesi [(| suggest an interpretation of 
their off-diagonal geometric phases in terms of Berry 
phases for single closed paths. In the second order case, 
these paths consist of the segments Ci, Gf-i, and Gik, 
where Gki geodesically connects the final point of Ck with 
the starting point of Cj, and vice versa for Gik (see fig. 
1 of ref. [6J). In the general non-Abelian case, how- 
ever, this interpretation is difficult to maintain. Apart 
from the special case when n\ = = ■ ■ ■ = n K , it 
is not possible to join the curves Cj lf . . . ,Cj„, due to 
the mismatch of dimensions, and thus the closure using 
geodesies is not applicable. This observation shows that 
the off-diagonal holonomies in general cannot be inter- 
preted as standard Wilczek-Zee quantum holonomies for 
closed paths [l4j, and that the off-diagonal holonomies 
therefore are genuinely new concepts associated with the 
evolution of quantum systems. Another consequence of 
the fact that we can have different n\ is that the rank of 

1 •" " cannot be larger than the smallest n; (see 2.17.8 
of ref. [29(), and thus may be less than n^. 

The non-Abelian character of the off-diagonal 
holonomies [C\, . . . , C K ] implies that they are not 
invariant under cyclic permutations of the indexes 
(Zi, . . . , Z K ). It may even be the case that two off-diagonal 
holonomies that differ only by a cyclic permutation have 
different rank, since the smallest ni only provides an up- 
per bound for the rank of 'y' 1 Furthermore, as is ex- 
emplified below it is also possible that -y' 1 ---'" may have 
path-dependent nodal points if K > 2. 



III. WHERE DID THE PHASE INFORMATION 

GO? 

As noted above, all strictly off-diagonal holonomies are 
undefined for cyclic evolutions, in analogy with the stan- 
dard cyclic off-diagonal geometric phases. Conversely, 
the noncyclic geometric phases are undefined when the 
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states at the initial and final points of the curves are or- 
thogonal, i.e., there exist nodal points where these phases 
are undefined. The hope to recover this lost phase infor- 
mation appears to have been one of the primary reasons 
for Manini and Pistolesi to introduce off-diagonal geo- 
metric phases. However, the issue concerning the pos- 
sible recovery of phase information was never explicitly 
investigated in ref. 6]. In the following we elucidate 
some aspects of this question in the case of non-Abelian 
off-diagonal holonomics. 

Let us first analyze what happens if the rank of some 
of the overlap matrices (!Fq\!F[) is greater than zero but 
less than their subspace dimension ni . This is a situation 
where the corresponding holonomies become partial [HI ; 
a phenomenon that has no counterpart in the Abelian 
case. To aid us in this analysis we introduce the unitary 
N x N matrix 



I <J . . . u 



(11) 



1)1 )))) 



It follows from unitarity that R(Stot) — N, where R(X) 
denotes the rank of matrix X. Furthermore, for every 
I = 1,...,J7, it holds that Efc erWt °' W = E*V fe tr ;fct = 
lnixnjj where l ni xm denotes the ni x ni identity ma- 
trix. This entails that (see 2.17.2 and 2.17.5 of ref. [H) 
£*-R(<* W ) > ni and Yj l R(o lh ) > n { . So, if i2(<r") = 
ni — n, then J2k^i R{") > n and J2k^i ^( CT ' fe ) > n - I n 
other words, when the overlap matrix (JqI^i) decreases 
by n in rank, the lower bound for the sum of the ranks 
of the matrices cr kl increases by the same amount. Thus, 
the "holonomy information" that is lost when the holon- 
omy of the curve Ci becomes partial is transferred to the 
matrices <r kl . 

A perhaps more significant question is whether the 
rank of the matrices -y' 1 ---'' 4 depends on the rank of the 
overlap matrix (!Fq\!F[) m a manner similar to what was 
discussed above for the matrices cr kl . One can demon- 
strate with a simple counterexample that no such relation 
exists. Assume r\ = 3 and n\ = ni = nj, — 2. Further- 
more, assume that <r 13 , <r 21 , and <x 32 are the 2x2 zero 
matrix, and 



lit 22 t 



_33t 



„12 ^23 „31 
(7 = (J (7 = 
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One may verify that the corresponding matrix Stot is 
unitary. In this example -y kl = 0, Vfc, I, and j klm = 
0, Vfc, I, m. Hence, although none of the overlap ma- 
trices (!Fq\!F[) are °1 ftd! rank, all strictly off-diagonal 
holonomies vanish, even those of higher order. 

Finally, we ask what happens if the matrices a 11 are 
zero for all I = l,...,r). We prove by reductio ad 
absurdum that at least one of the strictly off-diagonal 
holonomies must have nonzero rank. Assume cr 11 = 0, 
for I = 1,...,T), and 7 «*-«« = 0, for (h,...,l K ) G 1%, 
k > 2. Consider an arbitrary string (l±, ...,l v ) con- 
sisting of the integers 1 to 77, and with v > 2. If 



(h, . . . , l u ) G 12 {y < ??), then by assumption 7 Jl — = 
0. If (h,.-.,lv) ^ 1^) then take one of the small- 
est subsequences (l a , l a +\, ■ ■ ■ , lp) that begins and 
ends with the same number (i.e., lp = l a ). In this 
subsequence there is no other repetitions (otherwise 
there exists a smaller subsequence). It follows that 

over, (l a ,...,lp-i) G Ijg_ a _i. If - a = 1, then 
y«-'/3-i = cr lala = 0, otherwise <y l »-h-i = Q by 
assumption. We proceed by noting that Tr(S£ ot ) = 
j Tiry' 1- "'" = 0, for all v = 1,2,..., as a conse- 
quence of our assumptions. However, this cannot be the 
case since Stot is a unitary matrix. Therefore, our as- 
sumptions must be wrong, and at least one of the strictly 
off-diagonal holonomies must have nonzero rank. By the 
above findings we can conclude that the holonomy infor- 
mation lost in the the nodal points of the non-Abelian 
noncyclic holonomies indeed can be retained in some 
sense, but that the structure is much more intricate and 
rich than in the Abelian case, due to the existence of 
partial holonomies in the non-Abelian setting. 



IV. EXAMPLE 



We illustrate the off-diagonal holonomies by an ex- 
ample in the adiabatic context. We let the system 
evolve under the action of a slowly var ying Hamilto- 
nian. Let us consider the tripod system jl8l [30| mod- 
eled by the parameter-dependent four-state Hamiltonian 
H(s) — uj\e) (sin 9(s) cos (p(s)(0\ + sin#(s) siny?(s)(l| + 
cos 0(s)(a\) + h.c, exhibiting two nondegenerate 'bright' 
states l-Fjk) = {\B ± (s))} with energy ±w and a dou- 
bly degenerate 'dark' zero energy eigenspace \J-f) — 
{\D 1 (s)), \D 2 (s))}. Explicitly, we may choose 



\Di) 
\D 2 ) 



—=(\e) ± sin 9 cos ip\0) ± sin 9 sin ip\ 1) 
V2 V 



±cos9\a) J , 

cos 9 cos ip\Q) + cos 9 sin (p\l) — sin#|a), 
— sin^jO) + cos</?|l). 



12) 



Consider paths (0,0) — > (61, (pi) in parameter space 
(0,(p). For each such path the energy eigenstates de- 
fine paths C± in Q(i;l) and Cd in 0(4; 2). We obtain 
the geometric phase factors [C±] = 1 for 9\ ^ ir and 
U {2) [C ± ,C T ] = 1 for °i ^ 0. U {1) [C±] are undefined at 
0i = 7T and similarly U {2) [C±, C T ] at 9 X = 0. In ref. 0] it 
was shown that [Cd] is fully defined, except when the 
path ends at 9\ — 7r/2, where the holonomy becomes par- 
tial. The strictly off-diagonal holonomies (k = 2,3) in- 
volving the dark subspace are undefined when sin#i = 0. 
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U -p(o)-. ■ ■ -p(0)- y -f)(0) 



FIG. 1: Interferometric approach to obtain the holonomy 
Ug [Cjj, . . . ,C; K ]. The horizontal and vertical paths of the 
Mach-Zehnder interferometer correspond to the states |0) and 
1 1} , respectively. The particle have an internal degree of free- 
dom, on which the unitary operation U is applied. U can 
be generated either by adiabatic or nonadiabatic evolution, 
or approximated by filtering measurements. The particle en- 
ters the interferometer in path and in the internal state 
Ply (0)/niy . In path 1 a variable unitary operator V such that 
[V,P;(0)] = 0, V/, is applied, while in path we alternat- 
ingly apply U and a filtering corresponding to the projector 
|0> <0| <g> P ifc (0) + |l)(lj <g> 1, for k = 2, . . . k. A final U is ap- 
plied, followed by the filtering |0)(0| ®Pj x (0) + 11)(1| ® 1. The 
probability p to find the particle in path is measured after 
a second beam splitter. The unitary operator V is varied so 
as to maximize p. 



For sin 0i ^= 0, let Z = L cos[9(s)]tp(s)ds and we obtain 
U {2) [C ± ,C d ] = -C/ (3) [C±,C T ,C d ] 
= -UW[C±,C d ,C T ] -- 



cos((pi — Z) 
| cos((^i - Z)\ 

U {2) [C d ,C ± ] = -£/ (3) [C rf ,C±,C T ] 

I cos Z cos <pi sin Z cos (pi 



\ cos Z sin if i sinZsin^i 



(13) 



While U {2) [C d ,C±] and U i3) [C d ,C ± ,C T } are nonzero par- 
tial isometries, there are path-dependent nodal points 
of t/ (2) [C±,C d ], U^[C ± ,C T ,C d ], and U^[C ± ,C d ,C T ], 
namely where cos(</?i — Z) = 0. 



PHYSICAL REALIZATIONS 



Let us now examine some possible physical realiza- 
tions of U^[Ciy,. . . ,CjJ. Consider the Mach-Zehnder 
interferometer in fig. [T] with the two path states rep- 
resented by |0) and |1). We let the internal state of 
the particle (e.g., spin) be represented by the Hilbert 
space H = Hi(s) ffi • • • ffi H v (s), s £ [0, 1]. The total 
system is prepared in the state |0)(0| <g> Piy(0)/niy. We 
first apply a beam-splitter, followed by the unitary oper- 
ations |0)(0| + ®V and 1 0> <0 1 ®U+ |1)(1| ® 1, 
where V is a variable unitary operator assumed to be 
chosen such that [V, P/(0)] = for all I. Next, we per- 
form a filtering corresponding to the projection operator 



|0> <0j ®Pi 2 (0) + |1)(1| ® i, i.e., the particle is "removed" 
if it is found in path with its internal state outside 
subspace Ti.i 2 . Thereafter, we again apply the operator 
|0)(0|®C7 + |l)(l|(g)i, and the filtering |0)(0| <g> P* 3 (0) + 
|1)(1| (g) 1. This procedure is repeated until we have ap- 
plied the operator |0)(0| ® U + |1)(1| ® i, k times. After 
this, we apply a final filtering |0)(0| ® P h (0) + |1)(1| ® i, 
and recombine the two paths with a beam splitter. We 
finally measure the probability p to find the particle in 
path 0. 

First, we assume that the unitary operator U acting 
on Tt is caused by an adiabatic evolution of a time- 
dependent Hamiltonian with cigenspaces {Hz(s)}| ? _ 1 . 
This allows us to write U = J2i e where 4>i is the dy- 
namical phase 4>i = Jq 1 Ei(s)ds, and Ei(s) the eigenvalue 
corresponding to eigenspace TLi{s) of the Hamiltonian. 
The corresponding detection probability becomes 



P = 



1 11, , 

1 Tr(7;, ; 7, 

4 4fi,, Wii-U^. 



.) 



_IJLRe[e < £2=i**Tt(7 ,1 "- , "V t )l, 
2 n Zl 



where V lj = (1{{0)\V\1{ (0)). Note that V is a uni- 
tary matrix since [V, Pi ± (0)] = 0. By varying V we 
obtain the maximal detection probability when V = 
e »E£=i 4>kU^[Ciy,. . . ,CiJ. Hence, up to the dynamical 
phases we have found the holonomy. 

In the adiabatic setting there is in the general case no 
easy way to eliminate the dynamical phases. To avoid 
these problems we consider two alternative approaches 
to generate the unitary operator U. One alternative is 
to base the evolution entirely on filtering, where we ap- 
proximate the evolution in the spirit of ref . [22| . We be- 
gin with the same initial state, beam-splitter, and vari- 
able unitary V, as in the previous case. Next we ap- 
ply a sequence of filterings |0)(0| ® Ph{sj) + |1)(1| ® 1, 
where Sj form a discretization of the interval [0,1]. 
For the next step we apply the sequence of filterings 
1 0) <0 1 ® Pi 2 (sj) + |1)(1| (g) i, and we continue up l K . We 
finally apply |0)(0| ® P h (Q) + |1)(1| ® 1, followed by a 
beam splitter, and measure the probability to find the 
particle in path 0. One can show that the probability is 



(14) 



1 1 



1 

4 ' An h 
1 1 



Tr(7/ 1 ...i„7 i t 1 ..i re ) 
Re[Tr( 7 ' 1 -^y t )]. 



(15) 



Hence, as in eq. (|14[) . apart from the absence of dynam- 
ical phases. 

The second alternative that allows us to avoid the 
problem with dynamical phases is to use a nonadiabatic 
approach. Assume that the evolution is driven by the 
time-dependent Hamiltonian H(s), where now s is the 
time-parameter. We let the subspaces Tli(s) be evolving 
under H(s) according to the Schrodinger equation. In 
contrast to the adiabatic approach the subspaces TLi{s) 
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are in the general case not eigenspaces of H(s). We fur- 
thermore let {|J J '( s ))}i=i be smoothly parameterized or- 
thonormal bases of the subspaces Hi(s). We wish to find 
the unitary matrices Ui(s) such that the vectors 

ix?(»)> = £W»)>[tf«(«)]i* ( 16 ) 

i 

satisfy the Schrodinger equation id s \xf(s)) = 
H(s)\xi(s)) (ft — 1) with initial conditions 
|xf(0)) = \l k (0)). If we substitute eq. (I6J) into 
the Schrodinger equation we find that Ui(s) has 
to satisfy id s Ui(s) = iA t (s)Ui(s) + Ki(s)Ui(s), 
where [Ai(s)]fj = (d s P (s)\P(s)) and [Ki(s)]j>j — 
(P (s)\H(s)\P (s)) contain the geometrical and dynami- 
cal contributions, respectively, as was discussed in ref. 
[2lj . In order to get rid of the dynamical contribution 
without affecting the evolution of the subspaces we 
introduce the modified time-dependent Hamiltonian 

H(s)=H(s)-'£P l ( S )H(s)P l (s). (17) 
l=i 

The evolution of the subspaces Hi(s) are not affected by 
this modification since [H(s), Pi(s)] = [H(s), Pi{s)], VZ. 
We now wish to find the unitary matrices Ui(s) such that 
the vectors |xf(s)) = J2j W ( s )) [Ui(s)]jk satisfy the mod- 
ified Schrodinger equation idt\xi(s)) = H(s)\x k (s)) = 
[H(s) - Pi(s)H(s)Pi(s)]\xf(s)) with initial conditions 
|xf(0)) = \l k {0)). In this case we obtain d s Ui(s) = 
iAi(s)Ui(s). Hence, the solution Ui(s) = p e Io A '( s ') ds ' 
only depends on the geometric contribution 32] . The 
time-dependent Hamiltonian H(s) generates a unitary 
mapping from the initial state to the state at time s = 1 
that is given by £7(1) = J2i where 

r* = X3lxf(i)><xf(o)| 

k 

= ^[P e /o^] jfe |F(l))(/ fc (0)|. (18) 

jk 



This means that if we let £7 = £7(1) in the alternating 
procedure described above (see fig. [T]), then the proba- 
bility to detect the particle in path becomes as in eq. 
(p~5|) . and is maximized when V = U^[Ci lt . . . , C;J. 



VI. CONCLUSION 



Noncyclic evolution of quantum systems may lead to 
well-defined off-diagonal holonomies that involve more 
than one subspace of Hilbert space. These holonomies 
reduce to the off-diagonal geometric phases in ref. @ for 
one-dimensional subspaces. The off-diagonal holonomies 
are undefined for cyclic evolution but must contain mem- 
bers of nonzero rank when all the standard holonomies 
are undefined. While the nodal point structure of the 
holonomy for an open continuous path [24| can only de- 
pend on the end-points of the path, this structure can 
be path-dependent in the off-diagonal case. Further- 
more, we have put forward physical realizations of the 
off-diagonal holonomies in the context of adiabatic evolu- 
tion and interferometry that may open up the possibility 
to test these quantities experimentally. 
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